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ABSTRACT 


Torsional shafts made from fibrous composite materials 
are designed for minimum weight. Multiple loading conditions 
are considered. Design variables are the shaft inner diameter, 
the ply thicknesses, and the ply orientations. Design con- 
Straints include limits on ply strain, displacement, frequency, 
and Euler and shell buckling. The design task is solved as a 
numerical optimization problem. Examples are presented to 


demonstrate the method. 
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Lt. UNSROOUCTION 


Optimization and composite materials are two fields that 
have seen considerable developments in recent years. Advance- 
ments in these technologies have provided the design engineer 
with new dimensions, the depths of which are yet to be explored. 

The major material characteristics that are advantageously 
cultivated in composites are the high strength and stiffness 
Bomweieht ratios. The weight reduction is so significant that 
the aircraft and aerospace industries have concentrated major 
research efforts in this field. Investigations in improving 
the resistance of composites to environmental factors such as 
temperature, corrosion, and wear are continuing and findings 
affirm the contention that the future of composites is bright. 
However, other industries are slow in utilizing these materials 
in basic machinery elements, mainly due to a major drawback; 
the high cost of manufacturing and processing. It is in this 
vein that the importance of an optimized design becomes apparent. 
Advancement in fabrication and processing of composite materials 
coupled with the capability to optimize designs may well lead 
composites into a more competitive market. 

The application of numerical optimization techniques to 
the design of composite shafts is demonstrated here. The effects 
of a small mass imbalance on the deflection of a composite 


Shaft with synchronous whirl are included. 





Section 2 defines the scope and limitations of the analysis 
amd Optimization. 
Section 3 describes the analysis of composite shafts, 
mm@ecluding: 
A. Composite properties 
Be Shaft section properties 
fon sCommostce Latlure Criteria, ply strain calculations 
D. Displacement and vibration analysis 
E. Summary of failure modes considered 
Section 4 presents the optimization results [Refs. l, 2]. 
Three design examples are shown. 
Section 5 contains the conclusions drawn from this study 
and recommendations for future investigations. 
Appendix A contains a description of the analysis portion 


of the FORTRAN program used for composite drive shaft design. 





fil. S@@PE AND LIMITATIONS 


The analysis presented here applies to the shafts made of 
isotropic or composite materials with pinned-pinned end con- 
ditions. Basic shaft design formulas which include the effects 
of a small mass imbalance are used. The shaft is assumed to 
rotate with synchronous whirl. 

iiMewOeLamization iS constrained to subcritical speeds only. 
Supercritical speeds introduce feasible design regions that 
are disjoint from the primary design space. [In this investi- 
gation, these secondary areas are considered infeasible. 

Since the optimization tends to a large radius-small thick- 
ness design, radial stresses are neglected. However, buckling 
as a thin cylinder due to torsion or compression is considered. 
The equations used for these failure criteria do not include 
dymamic effects. A constant torque, centrifugal force, and 
axial force are considered in the formulas [Ref. 3]. 

The analysis incorporates the ability to design a shaft 
that can be used in two or more loading conditions. lifes 
example, a shaft may be designed to transmit 50 HP at 300 RPM 
Mee an axial load of 2000 lbs. as well as to transmit 150 HP 


MoO RPM with a LOOO Ibs. axial load. 


dels 





TEE. “Agee ysrs 


The problem considered here is a shaft required to transmit 


a specified horsepower at a given speed. Basic equations are 


used throughout the analysis. 


The inch-pound-second system of units (IPS) is used in the 


computer program analysis. However, by changing the equation 


for torque and the acceleration due to gravity to reflect the 


International System of Units, the program may be used with 


Pemsistent SI units. 


The assumptions used in the analysis are: 


loe 


+~ 


The plys are made of orthotropic material, but need 

not all be the same material. 

oG@euc ture: 

a. Plane stress in radial direction 

b. Plys are at different orientations and have 
different thicknesses. 

c. Cross-section is constant over entire length 


Loads: (see Figure 1) 


aoe eax al sogice 
IB) e Constant torque 
©. @G@entrifucal force 


Failure modes: 


Typical ply orientation and composite layup are 


Sijewn in Figure 2, 
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Stress element (Figure 2) 





F T 
dW 
Center of mass 
ply #n 
ply #2 
ply #1 


Elastic axis 


Peeters  wpplted Loads: aAxial load (F), 
Mass imbalance (e), Torque (T). 
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Composite layup 


ewe ee imeecal ply orientation and composite -avup. 
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a. SMaerinum ply strain limits 
b. Maximum displacement 
ec. First fundamental frequency 
d. Shell bwekling (axial and torsional) 
e. Euler buckling (contained within the displace- 
ment constraint) 
Ze COMPOSLTE PROPERTIES 
Stress is a measure of internal force within a body. This 
together with elastic constants are the variables for the 
determination of stiffness and strength of a material. The 
mechanism of deformation and failure are interpreted in terms 
of the state of stress and strain. When discussing stress it 
is usually considered to be the average stress over some physi- 
cal dimension. 
In Composites, three levels of average stress are considered: 
ga eromecmamtcal or local stress is that calculation 
based on distinct, continuous phases of fiber, 
matrix and, in some cases, the interphase and voids. 
b. Ply stress is that calculation based on assumed 
homogeneity within each ply or ply group where the 
fiber and matrix are smeared and no longer recog- 
nized as distinct phases. 
c. Laminate stress resultant N or moment M is based on 
an average of ply stresses across the thickness of 


a laminate. The individual plies are smeared. 


t- 
1 





ivemeeieettlctemes, Or material constants, of the stress- 
Strain relations can be packaged in set of engineering con- 
stants, compliance components or modulus components. 

a. MODULUS is used to calculate the stress from strain. 
This is the basic set needed for the stiffness of 
multidirectional laminates. 

Be COMPLTIANGE is used to calculate the strain from 
Stress. This is the set needed for the calculation 
of engineering constants. 

c. ENGINEERING CONSTANTS are the carryover from the 
conventional materials. Designers often feel more 
comfortable working with the familiar engineering 
constants. 

d. There is a direct relationship between the modulus 
and compliance. One is the inverse of the other. 

The stiffness of unidirectional composites is governed by 
the same stress-strain relation that is valid for isotropic 
materials. However the number of independent consitants are 
four for planar composites as compared to two for isotropic 
materials. The stiffness of isotropic materials can be repre- 
sented by the Young's modulus alone because the Poisson's ratio 
is the same all directions. 

For composite materials, Poisson's ratios are not bounded 
and can have very significant effect. Young's modulus alone 
1s not sufficient to describe the stiffness of composite 


Poe ees cise not SUtELcCient for isotropic materials either. 
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Pieeory sStitiness related to the global coordinate system 


are {Ref. 4]: 


8 4 22 a 2 
“eee ta Oe Oh TQ, + 4m n Q.. 


“ag nx yy 


_ 4 4 Zoe ae 
=n Q + m Q + 2m i Se + 4mon OP 


R59 ox yy 


z 2 Dee 


Dak Dae) 
=— — 
Qi mn Qx + mio Shows + (m +n ) 4Gmon Ca 





lage 2 Ae 2. oe 
= — + = 
266 mon Oe + min Oly 2m on , (m-<-n ) aes 
3 3 ee, Jee 
= —_ —_ — =< 
Q16 m nQ.x mn oe, + (mn -m > Cie 2(mn -m n)Q.. 
3 3 3 3 3 S| 
— — —_ + = 
Q5¢ mn Ore m nQy + (mn ~-mn Jee 2(m~> n-mn a 
where 
m= cos @ 
n = sin 9 
eg 
oe r i-v = 
= 
Q eae 
yy 1 oF 
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ireis Z te 
and Que? Quy? ae Olas On-axis modulus components referred to 
the local coordinate system 
Q. 53 Uti-a<ais Modulus components; i,J = 1,2,6 
MeeeonArFT SECTION PROPERTIES 
Axial stiffness, bending stiffness, and torsional stiffness 
are smeared properties. 
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He “Araraiwrst rt inesS (AE) 


Ninny 
AE = (Q Ge.) 
“el fae Ke k 
Des Bending stiffness (EI) 
NPLY 
Core), 1 C32) 
wey ok k 
eo, bhosstonralbustafttness. (GJ) 
NELLY 
GJ = > (Qi ute Cont) 
k=1 
where 
Qa the component of modulus my oy 0 
TL : the ply moment of inertia a inal esl i al Bop 
Jy chem piy polar Moment Of Inertia ; k ="1,NPLY 
AL. ;: the ply cross-sectional area oar 1 NPEY 
NPLY : the number of composite plys. 


meee COMPOSITE FAILURE CRITERIA 

For the determination of strength of any material it Is 
usual practice to estimate the stress at the time and location 
when failure occurs. In the case of conventional materials it 
is needed only to determine the maximum tensile, compressive, 
or shear stress or the principal stresses and then some obser- 
vation can be made avout the failure mechanism. This process 


is relatively straightforward because isotropic materials have 





Mompreterenmtial orientation and usually one strength constant 
will suffice. The isotropic material is essentially a one- 
Memensional or ome constant material. The Young's modulus for 
stiffness will suffice because Poisson's ratio is the same in 
all directions, and uniaxial tensile strength will also suffice 
because the shear strength is taken to be about 50 to 60 percent 
of the tensile strength. 

For composite materials, however, the on-constant approach 
for stiffness or for strength is no longer adequate. Four 
elastic constants are needed for the stiffness of a planer ele- 
Memc. Six constants for the strength of unidirectional composites 
are needed. Unidirectional composites have directionally depen- 
dent Strengths. The longitudinal strength can be twenty 
times that of the transverse and shear strengths; so for any 
State of stress, all three stress components must be examined 
before a judgement on the cause of failure can be made. The 
Specific stress component that is responsible for the failure 
may not be easily identified. Probably all three components 
are responsible. The effect of combined stresses must be sys- 
tematically determined and can be regarded as a way of life 
for composites. 

For composite materials, there is needed a failure criter- 
lan for the unidirectional plies. The strength of a laminated 
composite will be based on the strength of the individual plies 
within a laminate. Successive ply failures are expected as the 
@meetea LOad to a laminate increase. The first ply failure 


mmerotubowea by Gther ply failures until the last ply failure 


9 





which would be the ultimate failure of the laminate. The ply 
Seress and ply strain* calculations arewintended for strength 
determination. There are two popular approaches for failure 
criteria of unidirectional composites. They are each based on 
mem on—axis stress or strain as the basic variable [Ref. 4]. 


Jigs ine Maxsamum Stress Criteria 


Ox S ox 
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Failure occurs when one of the qualities is met. 


Ze. ie miascimum Strain Criteria 
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Failure occurs when one of the equalities is met. 
The maximum stress and strain criteria are not the same. 
Only when Poisson’s ratio of the unidirectional materials is 
mero, the criteria become identical. Conceptually they are 
Similar. Each component of stress or strain has its own cri- 
terion and is not affected by the other components. There is 


assumed to be no interaction. 


mech @em ctidy the ply strain failure criteria was used. 


D. DISPLACEMENT AND VIBRATION ANALYSIS 


In order to do the displacement and vibration analysis the 


deflected shape of beam is assumed as: 


20 





i of sin A (see Figure 3) 


where the beam is assumed to be simply supported. The deflec- 
Men is maximum at x = 1/2 [Ref. 3] 
4 K.e + K 
_ ye: 3 ¢ (ae eee 
a 384 (3D.1) 


2.4 22 


where: 
Ry = paw’ 
Ko face 
f = mass per unit volume 
A = cross sectional area 
® = shaft speed is radians per second 
g = acceleration due to gravity 
e = eccentricity of mass with respect to the axis of 
ECtEAeToOn 
2 = shaft length 
EI = smeared bending stiffness 


F = axial load 
Soest chamactLeristics of this equation are of particular 
mee rest. If there is no rotation or axial load, the formula 
reduces to the classical equation for maximum beam deflection 


under its own weight. 





4 
; 2K4 9 


max ~ 384ET eo 


Mewthere is no axial load, the equation becomes unstable 
4 
(deflection goes infinity) when EI = K, (4) OER sobvene Lor 


Ene critical speed, 


(3DR) 





2 EL 
pa 


which defines the first fundamental frequency of a simply 
Supported shaft. 


If there is no rotation, instabilty arises when 
OR a oboe) 
Beeognizable as Euler's column buckling criteria. 


Meee SOLn rotation and axial load, the critical speed 


becomes 


(2ipe ay) 





To avoid the instability regions, the denominator of 
equation (3D.1) must be given a lower bound greater than zero 
Meommt A in Fisure 4). Additionally, a maximum deflection has 


/ 


to be imposed to prevent computer overflow (point B in Figure 4). 


a2 





These limits have no relevance with the deflection constraint 
Mero in the design optimization process (point C in Figure 4), 
but are used only to prevent numerical ill-conditioning on the 


computer. 


Bee oUMMARY OF FAILURE MODES CONSIDERED 

Mabbure Cr2lteria serve important runctions in the design 
mmemesizing of composite laminates. They should provide a con- 
venient framework or model for mathematical operations. The 
framework should remain the same for different definitions of 
Mmpbures, such as the ultimate strength, the proportional limit, 
yielding, endurance limit, or a working stress based on design 
Or reliability considerations. Failures in composite materials 
involve many modes, such as fiber failures, matrix failures, 
interfacial failures, delamination, and buckling. Furthermore, 
the various modes interact and can occur concurrently or 
sequentially. 

MECCOCNPLY)+4) modes of failure are considered in this 
investigation. The 6(NPLY) (NLC) give the maximum strain limits 
gm Longitudinal compressive strain, longitudinal tensile strain, 
transverse compressive strain, transverse tensile strain, nega- 
tive shear strain, and positive shear strain. Failure is assumed 
memteply strain multiplied by a safety factor is greater than 
memescorresponding strain limit. 

Pe timit om the deflection is imposed. A common practice is 


PO Specify a maximum deflection in incnes per foot of shaft 


[2 
Lo 








Pinned-pinned end condition 


Bmore 5. Assumed deflected shape of the shaft. 


Deflection 
B \ 
C a a 
Ei 
A Denominator 


_ Figure 4. Limits onthe deflection. The solid curve 
mepeeoctics the deflection values used in 
the analysis. | 
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Memech. For example, as used in this report, 


df. = 0.005 inch/foot length 
max 


Banee the optimdzation may lead to a thin cylinder design, 
Baeckling due to torsion and compression are considered. The 
Metical torque and critical compressive stress for buckling of 


Meee cylinders are [Refs. 5, 6]: 


— _ ri¥2Q 
T = sara E rt ~ Baila: rt? 


Cc ate 3(1-y2)0°75 7 


Et atte: 
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Experimental data indicate that actual failure of the cylinders 
usually occurs below fifty percent of the values calculated 
from these equations. Therefore a large factor of safety is 
needed here. 

Lastly, failure is assumed if a specified percentage of 
the critical speed as calculated by equation (3D.5) is less than 
the shaft speed. 

The factors of safety used in this study are listed in 
the following: 


FALLURE MODE SAFETY FACTOR 
Strain 2 
Torsion Buckling 5 
Saecr!l Bwekling 5 


So) 


Frequency 


Deflection i 





Mm. OPDIMEZATION RESULTS 


There are a number of optimization programs available, 
each one using different technqiues in locating the desired 
optimum design. COPES/CONMIN [Refs. 1,2] is a versatile pro- 
gram that may be used for sensitivity analysis as well as an 
Optimization tool. Provided with a user supplied analysis pro- 
gram (Subroutine ANALIZ) where the objective function, con- 
straints and other relevant parameters are calculated, it 
determines a usable and feasible sector from which a search 
direction is chosen. This choice of search direction is made 
by using such information as the gradients of the objective 
Mamet ion and gradients of active and violated constraints. This 
iterative process of minimizing/maximizing the objective func- 
tion by changing the design variables is terminated when no 
further improvement can be made. 

Ma this study, the objeetive function to be minimized is 
the weight, with the shaft inner diameter andthe ply thicknesses 
and orientations as design variables. The analysis subroutine 
incorporates the analysis of composite shafts described in the 
previous sections. Several design examples are shown to demon- 
Strate the effects of the different parameters on the optimized 
design. For the examples considered here, the loads are listed 
in Table 1 together with the optimum weight and the constraints 


are listed in Table 2. 





Table 1: Optimization Examples 






Optimum 
Weight 

(lbs) 
OF 









Material 


G/E eS 
Steel ZO. 5 
J 
G/E 10.8 | 
| 
Steel 24822.0 
G/E S72. 5 


Examples No. 1 and No. 2 are designs for one load condition. 
Shaft lengths are 120 inches here. 

Examples No. 3 and No. 4 are designs for two load conditions. 
Shaft lengths are 60 inches here. 

Examples No. 5 and No. 6 are designs for two load conditions. 
Shaft lengths are 240 inches here. 


The results of optimization are presented in Figures 5-16. 


to 
ae | 





Table 2: Design Constraints 


Number Gone trait 


For NLC = 1 


Meenrough 6(NPLY) Maximum Strain Limits 
6(NPLY)+1 Buckling due tn Torque 

Oren e LY )+2 Buckling due to Compression 
eReNPLY)+3 Speed 

6(NPLY)+4 Deflection 


For NLC = 2 


eeereuY)+5 through 1L2(NPLY)+4 Maximum Strain Limits 
meaNPLY) +5 Buckling due to Torque 

me NPLY)+6 Buckling due to Compression 
L2(NPLY) +7 Speed 

meCNPLY)+8 Deflection 


Rte or NLC >? 2 
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COMPOSITE CRIVESRAFT OUTPUT 
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Figure 6. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


With all other parameters held constant, an increase in 
speed drives the design to a large inside diameter while 
reducing the thickness considerably. This trend is so strong 
that a reasonable lower bound on the thickness must be given 
to prevent numerical ill-conditioning. 

meme introduetion of eccentricity lets the optimization 
produce larger values for the design variables. 

The effect of axial load is reflected on the deflection. 

If very small deflections are allowed, the axial load effects 
are reduced. 

The versatility of the program should not be lost in the 
emphasis on composite shaft applications. The objective func- 
tion need not be confined to the weight. If the outer dimen- 
Sion is critical to the design, one may set a specific value on 
the diameter and let all other parameters vary during optimization. 

In summary, numerical optimization provides an efficient 
and effective way of creating composite shaft designs. 

The study has shown the feasibility of using numerical 
Optimization techniques in the designof composite shafts within 
the limitations imposed in the analysis. Further studies on 
the same design field may be pursued by aiminating some of these 
limitations. For example, the unfeasible region defined bv 
the speed constraint extends only to some discrete distance 


and a Supercritical feasible design area does exist. It will 
AL 





be interesting to develop optimization techniques to accomodate 
this disjoint design field and to investigate the peculiari- 
ties this will introduce to the anlaysis. 

Another investigation may be concerned with shafts with 
variable thickness and diameter along its length. The inclu- 
Sion of radial stresses in the calculations, the use of finite 
element methods in three dimensional stress analysis of the 
peeblem, and vibration considerations are but some of the aspects 
of research that may be useful in future anlaysis of composite 


materials. 
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APPENDIX A 


COMPUTER CODE USAGE 


The FORTRAN program used in this investigation is described 
here. Subroutine ANALIZ provides the basic analysis used in 
the optimization. It reads the initial design descriptions and 
M@eevcutlates the values of the objective function, constraints, 
and all other parameters necessary to solve the analysis 
problem. COPES/CONMIN updates the design to minimize/maximize 
the objective function, iterating until no further improvement 
meeene objective function is possible without violating one of 
the constraints. This process is shown in the flow diagran, 
Figure A-l. 

Mie global catalog given in Table 3 lists the locations, 
FORTRAN name, mathematical symbol and description of the para- 
meters used in the optimization. These parameters are contained 
in the labeled COMMON block, GLOBCM, which is accessed by 
@eeno for optimization. 

There are seven data groups to describe the initial design. 
The first data group contains the title. One card is used here. 
The second data group contains the number of plies, NPLY, the 
number of different materials, NMT, and the number of load con- 
ditions, NLC. One card is used, here. The third data group 
contains the initial inside diameter, shaft length, magnitude 
memmrass imbalance, fraction of shait radius at which the mass 


Mmmmerbance acts, and the acceleration of gravity. One card is 


es 
+ + 





used here. The fourth data group describes the thicknesses 
and the orientations of plies and material type used in the 
analysis. NPLY cards are used here. The fifth data group 
describes the ply longitudinal modulus, EL, ply transverse 
modulus, ET, ply shear modulus, GLT, ply major Poisson's ratio, 
PRLT, and the specific weight, RHO. The ply transverse Poisson's 
ratio, PRTL, is calculated internally. NMT cards are used here. 
The sixth datagroup specifies the ply longitudinal compressive 
oeeeein limit, EPLC, ply longitudinal tensile strain limit, 
EPLT, ply transverse compressive strain limit, EPTC, ply trans- 
verse tensile strainlimit, EPTT, and ply maximum shear strain 
limit, GMLT. The compressive strain limits are negative num=- 
bers, the tensile and the maximum shear strain limits are 
positive numbers. NMT cards are used here. The seventh data 
group describes the torque, moment, axial load, horsepower, 
and shaft speed in RPM for each loading condition. NLC cards 
are used here. 

The horsepower, torque, and shaft speed in RPM are inter- 
dependent values. [If two of them are known, the third one is 
calculated internally. The following table summarizes the 


required information for Subroutine ANALIZ. 
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Data 
Group 


Table 3: Input Data 


lntormae ion 


Title 

NPLY, NMT, NLC 

Digit. “Abe wnO, EC; GRAV 
iC) eto), MIYP (1) 

Pe eeteeechl, PRLT, RHO 

Polo, Erb i ere, MPrIt, GMLT 


T, M, F, HP, RPM 


-~ 
i 


Format 


20A4 


S20 


pu 0. 


2 ee O. 


SFi0. 


ea 00) 


3) 035 Oe 


2 0 





Call 


SUBROUTINE ANALIZ (ICALC) 


Flow Diagram 


ICALC = 1 


Subroutine SHFTO6 

Read Initial Dimensions 
Read Material Properties 
Read Design Conditions 


Pra cecatet lope 


ICALC = 2 


Calculate Geometric Parameters 


Gave 


Cagle. 


a 


Cali 


Gaile 


Call 


Subrowrtine sire, 

Calculate Objective Function 
SUubmRomeine SHETO3 

Cameulatesshatt Section Properties 
Subroutine SHEFTO9 

Calculate All Constraints 
Subroutine SHFT10 

Calculate Critical Frequency 
Subroutine SHEFTLIL 


Calculate Maximum Deflection 


ICALC = 3 


Subroutine SHFTO8 


Print Analysis Results 


{~ 
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OPTIMIZATION FLOW DIAGRAM 


INPUT 


CORES _——. CONMIN 


SUBROUTINE 


ANALI2Z 





OUReUT 


Freese -\- | 





Global 


mocation 


i 

2 

3 

4-10 
iLalysags 0 
31-50 
Bel > > 1 


Table 4: Global Catalog of Parameters 


FORTRAN 


Name 


VOL 
DIAMI 


WGHT 


TIN (20) 


THND (20) 


G (500) 


4§ 


MATH 


Symbol 


V 


d 


Dern con 


Volume 

loner Daamneter 
Weight 

Dummy Storage 
Thickness 


Orientation 
(Degree) 


Constraints 
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